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ABSTRACT 


This  document  describes  the  development  of  Unitless  Mathematical 
Models  for  on-board  flight  software.  The  development  is  part  of  the  NRL 
Software  Cost  Reduction  Program.  A  subset  of  the  equations  existing  in  the 
operational  flight  program  were  examined  to  determine  existing  units  and  basic 
assumptions.  The  equations  were  then  re-derived  in  unitless  form  with  clearly 
stated  assumptions;  this  is  the  form  which  is  most  useful  to  modular  design. 
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SECTION  1 
INTRODUCTION 


1.1  Background 

The  concepts  of  Modular  Mathematics  and  Information  Hiding  are  now 
being  applied  to  software  for  the  A-7  aircraft  on-board  computer.  This 
demonstrat ion  program  has  been  active  for  several  years.  However,  to  fully 
apply  these  concepts,  it  has  become  evident  that  a  detailed  analysis  of  major 
portions  of  the  basic  mathematics  is  required.  There  are  several  reasons  for 
this.  First,  there  are  a  number  of  equations  within  the  Operational  Flight 
Program  (OFP)  that  contain  constants  with  units  which  are  unknown  or  do  not 
conform  to  the  basic  data  types.  It  is  not  evident  by  inspection  how  these 
constants  were  determined  or  on  what  they  might  depend.  That  is,  it  is  not 
apparent  how  such  constants  depend  on  basic  quantities  with  known  physical 
meanings.  Second,  the  mathematics  also  depend  on  certain  basic  assumptions 
not  readily  discernible  in  the  form  in  which  the  equations  are  presented.  The 
nature  of  these  assumptions  must  be  known  so  that  a  given  model  may  be  placed 
into  a  specific  module.  Then,  if  a  basic  assumption  changes,  the  software 
modification  should  be  confined  to  a  single  module.  Finally,  the  equations 
used  in  the  OFP  are  not  in  a  form  which  supports  modularity.  Instead,  the 
equations  produce  intermediate  products  which  do  not  conform  to  the 
established  basic  data  types. 

Modular  Mathematics  can  be  usefully  defined  as,  Mthe  statement  of 
the  governing  equations  in  terms  of  known  assumptions  using  'universal' 
constants  and  in  a  unitless  formulation."  Modular  Mathematics  supports  the 
principles  of  Information  Hiding  by  allowing  the  user  to  be  concerned  only 
with  constants  of  known  physical  meaning  and  not  with  the  actual  units  used. 
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Aside  from  being  a  necessary  condition  for  the  employment  of  the  principle 
of  Information  Hiding,  Modular  Mathematics  is  useful  in  its  own  right.  It 
provides  a  cost  reduction  during  the  development  phase  by  highlighting 
assumptions  inherent  in  the  mathematical  statement  of  the  problem.  Thus,  an 
algorithm  that  optimizes  the  time,  accuracy,  and  storage  requirements  can  be 
selected.  A  carefully  compiled  set  of  equations  which  are  rigorously  derived 
and  fully  documented  should  be  highly  reusable.  For  example,  a  well-written 
and  documented  subprogram  will  find  frequent  reuse  to  meet  similar  needs  in 
other  systems.  It  will  also  serve  as  a  model  and  point  of  departure  for 
improved  versions  of  the  same  subprogram  or  other  subprograms  developed  to 
meet  related  tasks.  Thus,  a  well-supported  mathematical  development  will 
serve  the  user  who  only  needs  the  mathematics  to  solve  a  similar  problem  or 
needs  a  good  example  to  follow  in  other  mathematical  developments,  in  perhaps 
unrelated  areas.  Once  a  set  of  software  is  developed,  the  modular  form  of  the 
mathematics  can  provide  substantial  aid  during  program  maintenance.  A  change 
in  device  or  system  specification  can  be  relatively  quickly  traced  to  the 
relevant  equations  and  from  there  to  the  point  requiring  modification  in  the 
derivation  or  change  in  the  constants.  Then  the  change  can  be  made  rapidly 
and  with  little  danger  of  any  widespread,  undesired  effect.  Substantial  cost 
savings  can  be  realized  by  not  having  to  reproduce  the  mathematics  from  first 
principles  for  each  change  in  device  or  system  specification. 

1.2  Techniques  for  Generating  Unitless  Models 

A  Unitless  Model  is  not  only  a  mathematical  prescription  which  is 
free  of  unknown  constants  or  operations  with  inappropriate  data  types;  it  is 
the  entire  derivation  behind  that  equation,  including  the  origin  of  the 
relevant  equation  and  the  sequence  of  mathematical  operations  that  follow  to 
arrive  at  the  final  useful  form.  The  trail  includes  all  assumptions  made  and 
all  constants  identified  and  introduced.  In  an  ideal  situation,  the  model 
could  be  constructed  from  scratch.  Consideration  could  be  given  to  such 
issues  as  accuracy,  machine  time  required,  storage  required,  error  propaga¬ 
tion,  and  simplicity  of  form.  In  many  cases,  the  mathematics  already  exists, 
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so  the  origin  and  derivation  must  be  determined  from  the  form  of  the  equation 
and  knowledge  of  the  relevant  processes. 


The  mathematical  analysis  Is  approached  by  considering  a  detailed 
statement  of  the  problem  taken  from  the  description  of  the  relevant  module. 
The  module  description  may  not  contain  enough  detailed  information  to  set  up 
entirely  the  mathematical  problem,  and  the  description  must,  therefore,  be 
expanded  to  include  a  sequence  of  steps  or  operations  that  will  accomplish 
this  task.  The  new  problem  statement  must  be  examined  to  determine  if  it 
still  adheres  to  the  basic  tenents  of  the  module  and  if  it  is  a  clear  and 
concise  statement.  In  many  cases,  several  derivations  are  possible.  The 
selection  between  them  may  be  performed  at  this  point  in  the  development  if 
there  is  a  clear,  compelling,  distinction  between  them.  If  the  decision  is 
not  made  here,  then  the  alternative  forms  may  be  carried  forward  for  a 
decision  at  a  later  point  in  the  development,  based  on  additional  information. 

The  first  step  is  to  select  a  starting  point  for  the  mathematical 
development.  It  should  be  as  general  and  free  from  restrictive  assumptions  as 
possible.  Consider  a  ballistics  problem,  for  example.  Rather  then  stating 
that  a  force  is  equal  to  mass  times  acceleration,  we  could  start  with  the  more 
general  case  of  conservation  of  momentum.  We  would  then  argue  that  since  mass 
does  not  change,  the  result,  F~ma,  is  applicable.  By  starting  from  the  more 
general  statement,  we  have  highlighted  an  assumption  (i.e.,  that  mass  does  not 
change)  which  is  implicit  to  the  problem.  Thus,  in  selecting  a  point  of 
origin,  one  should  use  the  most  general  statement  possible.  Of  course,  this 
principle  can  be  carried  to  extreme.  Thus,  that  point  of  origin  should  be 
selected  that  includes  the  full  range  of  different  assumptions  encountered, 
yet  takes  advantage  of  assumptions  that  are  never  expected  to  change. 

The  mathematical  development  should  be  as  mathemat ically  rigorous  as 
is  consistent  with  the  manner  in  which  the  equations  are  used.  The  assump¬ 
tions  and  approximat Ions  made  during  the  development  must  be  justifiable  and 
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consistent  with  usage  of  the  final  form.  These  assumptions  should  be 
carefully  noted  in  the  development.  In  the  case  where  there  is  no  preexisting 
final  form  of  the  equation,  considerat ion  must  be  given  to  required  accuracy, 
storage,  timing,  and  other  parameters  in  selecting  the  end  point  of  the 
development.  Accuracy  considerations  require  that  the  equations  derived 
include  all  corrections  which  have  a  measurable  effect  somewhere  in  the  module 
result.  However,  effects  several  orders  of  magnitude  down  from  the  accuracy 
level  of  the  given  quantity  should  be  ignored.  Since  no  computer  offers  truly 
unlimited  resources  of  time  and  storage,  consideration  must  be  given  to 
minimizing  the  use  of  both.  The  traditional  programming  technique  to  reduce 
computational  time  is  to  store  as  many  constants  as  possible;  on  the  other 
hand,  to  minimize  storage,  constants  are  typically  recalculated.  In  the 
context  of  Mathematical  Modularity,  either  approach  must  be  used  with 
discretion.  Constants  selected  for  storage  must  have  a  recognizable  physical 
meaning  so  that  they  can  be  available  to  the  correct  module.  Further, 
constants  should  have  a  dimensionality  consistent  with  the  basic  data  types. 

Once  the  final  equation  is  derived,  it  must  be  posed  in  a  form  which 
is  consistent  with  the  basic  data  types.  The  basic  data  types  are  selected  to 
meet  the  needs  of  the  problem  at  hand.  They  are  the  minimum  set  which  meet 
the  needs  of  the  program  while  not  taxing  the  system.  The  final  form  of  the 
equation  should  include  only  those  operations  which  produce  one  of  the  basic 
data  types  as  the  end  product. 

1.3  Report  Organization 

This  report  documents  the  results  of  the  initial  phase  of  the 
project  to  generate  Unitless  Models  for  the  A-7  flight  software.  The  report 
is  divided  into  chapters  which  parallel  the  division  of  the  program  into 
modules.  Each  chapter  is  self-contained  and  includes  lists  of  symbols  and 
references  pertaining  to  that  chapter;  more  chapters  can  be  added  as 
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additional  modules  are  analyzed.  This  will  also  allow  corrections  to  be  made 
to  a  single  chapter  without  the  necessity  of  changing  anything  in  any  other 
chapter.  The  notation  and  references  used  are  consistent  throughout  the  text 
and  are  also  listed  separately  in  appendixes. 
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SECTION  2 
AIR  DATA  COMPUTER 


2.1  General 

The  Air  Data  Computer  (ADC)  processes  inputs  from  the  Pitot-static 
system  and  the  total  temperature  probe  to  provide  outputs  that  represent 
barometric  altitude,  Mach  number,  and  true  airspeed. 

2.2  Altitude  Correction  for  Non-AIMS-probe  and  AIMS-probe  systems 

Background 


Page  ADC- 13  of  Reference  1  on  the  ADC  gives  two  equations, 
alt  alt  +  [560.  *{ cmindex(mi)  -  .2)]  (2-1) 


and 


alt  :=  alt  *  {  1.  -  cmindex(mi )**2 

*f  .02032  -  ( (7.5l48E-7)*alt) ] } 

+C0NVERT( ! system  vertical  velocity!).  (2-2) 

These  equations  correct  the  measured  value  of  the  barometric  altitude  for 
errors  associated  with  the  location  of  the  Pitot-static  port. 


There  are  two  Pitot-static  systems  available  for  the  A-7  aircraft. 
The  first  has  the  static  ports  along  the  side  of  the  aircraft  just  below  the 
canopy  rail.  The  second  has  its  static  ports  located  on  the  Pitot  tube,  which 
is  an  L-shaped  tube  protruding  from  the  side  of  the  aircraft  just  aft  of  the 
radome.  The  second  system  is  referred  to  as  an  AIMS  (or  in  some  cases  an 
L-probe)  system,  and  the  first  as  a  non-AIMS  system  (or  non-L-probe  system). 
Equation  2-1  provides  correction  on  the  non-AIMS-probe  aircraft,  while 
Equation  2-2  refers  to  AIMS  aircraft.  Several  quantities  need  to  be  defined. 

Static  Pressure  -  The  absolute  pressure  of  the  still  air  surrounding  the 
aircraft . 

Static  Defect  -  The  difference  between  the  pressure  at  various  places  on 
the  aircraft  skin  and  the  free-stream  pressure. 

Position  Error  -  The  error  in  the  static  pressure  at  the  static  port  due 
to  static  defect. 

The  pressure  on  the  aircraft  skin  may  be  either  slightly  lower  or 
higher  than  the  free-stream  pressure.  The  static  defect  at  a  particular 
location  depends  on  speed,  angle  of  attack,  and  altitude.  To  obtain  a  sample 
of  static  air  on  a  moving  aircraft,  a  hole  (static  port)  or  a  series  of  holes 
can  be  drilled  in  a  plate  on  the  side  of  the  fuselage  or  on  the  side  of  a 
rod-like  probe  (called  a  Pitot-static  tube)  extending  into  the  free-air 
stream.  The  location  of  the  static  port  is  selected  by  wind-tunnel  tests  and 
by  tests  at  numerous  locations  on  the  actual  aircraft.  The  static  port  is 
connected  to  a  pressure  transducer  by  means  of  tubing.  The  time  constant  of 
the  cavity  formed  by  the  tubing  and  the  transducer  is  a  function  of  the 
viscosity  of  the  air,  cavity  volume,  static  pressure,  and  length  and  diameter 
of  the  connecting  tubing. 

The  time  required  for  a  change  in  static  pressure  to  be  sensed  by 
the  Pitot-static  system  on  a  maneuvering  aircraft  introduces  further  errors. 
An  additional  term  must  be  defined. 
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Pressure  Lag  Error  -  The  time  delay  between  the  change  in  static 

pressure  and  the  change  in  signal  from  the 
pressure  transducer. 

Thus,  there  is  no  pressure  lag  error  in  straight,  level,  unaccelerated  flight. 
Definitions  of  Symbols 

We  define  the  following  variables  and  constants  of  a  general  nature: 

G  force  of  gravity  at  mean  sea  level 

32.17405  ft/sec2 
9.80665  m/sec2 

H  scale  height  of  the  atmosphere  27672.24  ft. 

hp  barometric  altitude  corrected  for  position  and  pressure  lag 

error 

hp  ADC-sensed  barometric  altitude 

Mn  true  Mach  number 

PQ  standard  pressure 

29.92  in.  Hg. 

101332.27  N/m2 
2116.2  lb/ft2 

Ps 


static  pressure 


r 


t 


» 


t 


* 


l 


Pt  total  pressure  (sum  of  dynamic  and  static  pressure) 

q  dynamic  pressure  (P^  -  Ps) 

R  gas  constant  for  dry  air 

1716.5  ft2 /sec2 -°R 
287.04  m2/sec2-°K 

T0  standard  temperature 

518,688  °R 
288.16  °K 

Vvv  vertical  velocity  or  rate  of  change  in  altitude 

a  standard  temperature  lapse  rate  of  the  atmosphere 

0.00356616  °F/ft 
0.0065  °C/m 

cxj  true  angle  of  attack 

y  ratio  of  specific  heats 

1.4  for  air 

t  time  constant  associated  with  the  static  pressure  system 

Derivation 


The  measured  barometric  altitude  must  be  corrected  for  position 
error  and  pressure  lag  error.  The  functional  forms  of  these  quantities  are 
determined  by  wind  tunnel  and  flight  testing.  In  the  preceding  section,  it 
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was  noted  that  position  error  is  a  function  of  speed,  angle  of  attack,  and 
altitude,  while  pressure  lag  error  is  a  function  of  the  time  constant  of  the 
cavity.  This  can  be  expressed  in  the  form, 


h 

P 


f  (M  ,  a_,  h  )  +  f  (t,  V  )  , 
p  ^  n  T  pJ  t  ^  v v 


(2-3) 


where  fp  denotes  the  functional  form  of  position  error  and  fT  that  of 
pressure  lag  error.  The  exact  nature  of  both  of  these  functions  is  determined 
by  fitting  curves  through  the  experimental  flight  test  data. 


Now  consider  the  pressure  altitude  only  in  terms  of  the  parameters 
that  form  the  position  error,  fp,  or 


Then  a  Taylor  series  expansion  of  hp  (Mn,  ay)  about  Mn  =  Mn^  and 
ay  =  ayQ  yields 


h  fM 


aT)  =  h  (M  ,  aT  ) 
p  v  n  T  } 
K  o  o 


+  [M  -  M  1 
v  n  n  J 
o 


+  0  (AM^,  Aa^.^) 


(2-4) 
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We  let 


h 

P 


and  see  that  hp  can  be  treated  as  a  constant. 


Comparing  Equation  2-4  with  Equation  2-1,  we  see  that  the  former 
will  reduce  to  the  latter  if  the  position  error  is  only  a  function  of  Mach 
number,  and  if  a  constant,  kj  ,  is  defined  as 


The  resulting  equation  is 


h 

P 


(2-5) 


The  constants,  kx  and  Mn^,  are  empirically  determined  from  flight  test 
data. 


If  a  similar  comparison  is  made  between  Equation  2-2  and  Equation 
2-4,  no  formal  resemblance  Is  noted.  However,  if  the  error  in  the  measurement 
of  true  altitude  due  to  position  error  is  assumed  to  be  a  function  of  dynamic 
pressure  (q)  and  Mach  number  (Mn)  rather  than  speed,  angle  of  attack,  and 
altitude,  then  a  Taylor  series  expansion  of  hp  (q,  Mn)  about  q  =  qQ 
Mn  =  Mn<)  yields 
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Vq’MJ  =  hp(V  Mn  )  +  h  -  (—  hp  MrJ]  M 

H  K  o  \3q  /  V  Mn 


+  (Mn  -  Mn  )  —  hp  (q>  Mn^l 

«  H  P  )%■  Un 


+  0  (Aq2,  AM^2)  . 


(2-6) 


We  let 


h 

P 


9 


and  see  that 


can  be  treated  as  a  constant. 


Equation  2-6  must  be  related  to  f light-measurable  quanitites. 
Consider  Bernoulli’s  Equation, 


(2-7) 


and  the  defining  equation  for  the  standard  atmosphere, 


(2-8) 


2-7 


Using  hp  to  approximate  hp  in  Equation  2-8,  and  substituting  for  Ps  in 
Equation  2-7,  yieids 


from  which  it  is  seen  that  we  may  consider  q  =q  (Mn).  Thus,  by  the  chain 
ruie, 


3  h  3h  dq 

_ p  _  _ P  __ 

3  M  3  q  dM 

n  n 


and  so  Equation  2-6  can  be  written  as 


h  «  h  + 
P  P 


(q  -  qQ) 


+  (M  -  M  ) 
v  n  n 

o 


(2-10) 


Taking  derivatives  in  Equation  2-9,  we  have 


(2-11) 


Since  q0  is  defined  to  be  a  reference  dynamic  pressure  at  the 
point  where  the  correction  is  zero,  it  foiiows  that  q0  =  0  at  Mn  =  0.  In 
view  of  this,  and  using  Equations  2-9  and  2-11,  we  can  write  Equation  2-10  as 


2-8 


or 


+  Y  1 


(2-12) 


i  oth 

Since  I  - & 

I  T 

o 

theorem  may  be 


<  1  and  M  2  <  1  for  all  reasonable  values,  the  binomial 

IV  2  /  n  I 

used  to  generate  the  sums 


(2-13) 


(2-14) 


and 


1 


(2-15) 


Actually  applying  Equations  2-13,  2-14,  and  2-15  to  Equation  2-12, 
and  neglecting  resultant  terms  containing  a  or  for  j  _>  4,  yields 


h  «  h 
P  P 


H  2 
n 


(*1  -  K2 


v 


(2-16) 


2-9 


where 


> 


f 


t 


3 

3q 


0,0 


and 


K2  =  -  Kr 

RT 


The  form  of  Equation  2-16  now  matcht:^  that  of  Equation  2-2.  The  constants, 
Ki  and  depend  on  the  empirical"  determined  value  of  the  derivative. 

This  derivative  -  or  more  likely.,  the  values  of  and  K2  -  are  determined 
by  flight  and  wind  tunnel  testing.  The  constant  K-|  is  unitless,  while  K2 
has  units  of  length.  Rewritir g  Equation  2-16  in  a  slightly  different  form 
gives 


f 


» 


I 


The  constant,  H,  introduced  in  Equation  2-17  is  the  scale  height  of  the 
atmosphere, 

T  R 


To  five  decimal  places,  the  constants  are 


=  .02032 
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l 


and 


» 


I 


t 


.02080  . 


The  difference  in  the  values,  and  *  *s  due  t0  *he  ^act  that  the 
original  constants  were  empirically  determined.  A  simplified  form  of  Equation 
2-17  is  possible  if  a  single  value  is  used  for  both  constants.  This  result  is 


* 


t 


« 


I 
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h  »  h 

P  P 

with 

=  .02032  . 

The  next  step  is  to  include  the  effects  of  pressure  lag  error.  Assumption  of 
a  linear  dependence  of  altitude  correction  on  vertical  velocity,  Vvv,  gives 

f  =  -K0  V  +  c  .  (2-19) 

T  3  vv 


1  - 


M  2 
n 


(2-18) 


Since  there  is  no  correction  in  level  flight,  the  constant,  c,  is  zero.  The 
constant,  K3,  has  dimensions  of  time  and  is  the  time  constant  of  the 
cavity.  The  non-AIHS  system  is  not  corrected  for  the  cavity  time  constant, 
while  the  AIMS-system  correction  becomes 


(2-20) 
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* 
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» 


where  and  are  defined  as  above  and 

K3  =  ,4375  sec. 

The  equations  for  the  position  and  pressure  lag  errors  are  dependent  on  the 
Pitot-static  system  used,  the  physical  dimensions  of  that  system,  and  the 
standard  atmosphere  assumed. 

The  only  remaining  correction  is  the  offset  applied  to  the  value  of  altitude 
reported  by  the  pressure  transducer.  The  OFP  uses  an  offset  of  -1024  feet  (a 
convenient  power  of  two),  while  the  assumed  offsets  for  the  non-AIMS-probe 
system  and  the  AIMS-probe  system  are  -1020  feet  and  -1056  feet,  respectively. 
The  final  form  of  the  equation  for  non-AIMS-probe  should  then  be 


I 


» 


h  =  h  +  K  fM  -  M  )  -  K  .  (2-21) 

P  P  1  V  n  no;  2 


The  final  form  of  the  equation  for  the  AIMS-probe  is 


The  values  and  physical  meanings  of  the  constants  appearing  in  Equation  2-22 
are  given  in  Table  2.1. 


> 


» 


> 
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Table  2.1  LOCAL  CONSTANTS 


Physical  Meaning 


Dependencies 


instanteous  rate  of  change  altimeter 

of  ratio  of  sensed  baro¬ 
metric  altitudes  and  Mach 
number 


Non-AIMS  offset 


OFP 


Mach  number  offset 


empirical  testing 


empirical  testing 


empirical  testing 


time  constant 


empirical  testing 


Value/Units 
560  feet 

1020  feet 

.2  unitless 

.0232  unitless 

.02080  unitless 

.4375  sec 


AIMS  offset 


OFP 


1056  feet 


2.3  Correction  to  Altitude  for  Non-Standard  Sea  Level  Pressure 


Background 

Page  ADC- 13  of  Reference  1  on  the  Air  Data  Computer  gives  the 

equation, 


Pl:  =  cfeet_dist  {alt  +  [cingh(slp  -  ingh(29.92) )]  (2-23) 

*[ 924.87  -  .00635856*alt] }  . 

This  equation  corrects  for  a  non-standard  sea  level  pressure  (i.e.  a  pressure 
other  than  29.92  inches  Hg  at  sea  level,  and  set  into  the  computer  by  the 
pilot) . 


There  are  several  important  terms  that  are  relevant  to  the 
measurement  of  altitude. 

Tape  line  altitude  -  The  actual  altitude  that  an  hypothetical  ruler  or 

tape  would  measure  from  mean  sea  level  to  the 

aircraft.  This  altitude  is  that  which  would  be 
given  by  a  perfect  altimeter. 

Geopotential  altitude  -  The  measure  of  the  potential  energy  of  a  unit 

mass  at  this  point  relative  to  mean  sea  level. 

The  geopotential  altitude  is  equal  to  the  tape 
line  altitude  only  if  the  force  of  gravity  is  not 
a  function  of  altitude.  Geopotential  altitude 
must  be  considered  in  all  pressure  or  density 

calcuiat ions. 


if 


f 


f 


! 


» 
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Pressure  Altitude  -  The  altitude  given  by  an  aneroid  altimeter.  It 

is  a  measure  of  the  difference  in  pressure 
between  a  standard  datum  (usually  mean  sea  level) 
and  the  measured  static  pressure  converted  to 
distance  by  assuming  a  standard  atmosphere. 

The  problem  encountered  here  is  that  the  pressure  altitude  furnished  to  the 
ADC  is  determined  using  the  standard  atmosphere,  or  29.92  inches  Hg.  This 
pressure  may  or  may  not  correspond  to  that  encountered  at  mean  sea  level.  The 
required  altitude  is  that  determined  using  the  actual  pressure  at  mean  sea 
level . 

Definitions  of  Symbols 

We  define  the  following  variables  and  constants  of  a  general  nature: 

G  acceleration  due  to  gravity  at  mean  sea  level 

32.17405  ft /sec2 
9.80665  m/sec2 

g(z)  acceleration  due  to  gravity,  a  function  of  altitude 

hp  geopotential  altitude  (see  Figure  2.1) 

H  scale  height  of  the  atmosphere  27672.24  ft. 

Hph  pressure  altitude  corrected  for  actual  pressure  at  sea  level 
( see  F igure  2. 1 ) 

P0  standard  pressure  (see  Figure  2.1) 


» 


29.92  In  Hg 
101332.27  N/m2 
2116.2  lb/ft2 
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PRESSURE 


Ps  measured  value  of  the  static  pressure 

R  gas  constant  for  dry  air 

1716.5  ft2/sec2-°R 
287.04  m2/sec2-°K 

T0  standard  temperature 

518.688  °R 
288.16  °K 

z  tape  line  altitude,  i.e.,  the  actual  distance  from 

mean  sea  level  to  the  aircraft 

a  standard  temperature  lapse  rate  of  the  atmosphere 

0.00356616  °F/ft 
0.0065  °C/m 

AP  difference  between  standard  pressure  and  actual  pressure 

at  sea  level 


Derivation 


To  arrive  at  Equation  2-23,  consider  a  Taylor  series  expansion  about 
the  altitude  (hp)  computed  using  the  standard  pressure  datum  (PQ),  which 
yields 

3  h 

H  =  h  +  AP  — E  +  0(AP2 ) .  (2-24) 


2-17 


> 


» 


t 


If  terms  of  the  order  of  AP2  and  higher  are  neglected,  then  the  remaining 
problem  is  to  determine  3hp/8P0.  The  standard  atmosphere  is  given  by  the 
equation, 


(2-25) 


1 


Since  the  static  pressure,  Ps,  is  not  known  to  the  module,  it  is  necessary 
to  calculate  it  based  on  the  pressure  altitude  and  other  known  sea  level 
quantities.  The  derivation  begins  with  the  definition  of  the  temperature 
lapse  rate,  a ,  or 

-  =  -a  •  (2-26) 

dh 

P 


Both  sides  of  Equation  2-26  are  multiplied  by  G/RT  to  give 

G  dT  Gd 

RT  dh  RT 

P 


(2-27) 


Geopotentlal  altitude  is  defined  as  a  measure  of  the  gravitat ional  potential 
energy  of  a  unit  mass  at  a  point  relative  to  mean  sea  level.  It  can  be 
defined  in  differential  form  by  the  equation, 

Gdhp  =  g(z)dz.  (2-28) 

Equations  2-27  and  2-28  along  with  the  chain  rule  are  used  to  give 

G  dT  dz  Cot 

RT  dz  dh  RT 

P 


I 


I 
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or 


(2-29) 


ft 


G_  £T  g(z)ot 

RT  dz  RT 

Next,  use  is  made  of  the  following  assumptions:  (1)  the  air  is  dry,  (2)  the 
atmosphere  is  a  perfect  diatomic  gas,  and  (3)  the  atmosphere  is  in 
equilibrium. 

The  atmosphere  then  obeys  the  perfect  gas  law, 

P  =  pRT,  (2-30) 

where  p  is  the  density  of  the  atmosphere. 

For  the  case  of  hydrostatic  equilibrium,  it  can  be  assumed  that 

dP  ,  .  , 

—  =  -  pg(z).  (2-31) 

dz 

In  view  of  Equation  2-31,  Equation  2-29  can  be  written  (again  using  the  chain 
rule)  as 

£_  dT  dP  _  _  g(z)q 
RT  dP  dz  RT 


I 


» 


Use  of  these  results  and  Equation  2-31  yields 

-  "  ( -pg(z>)  =  .  Siili  , 

RT  dP  RT 


G  dT  _  _1 
RaT  dP  P 
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Thus  we  have 


G  dT 
fa  T 


dP 

P 


(2-32) 


Now,  integration  of  the  defining  equation  for  the  temperature  iapse  rate 
(Equation  2-26)  gives 

T  h 

/  dT  =  -  a  /  P  dhn 


or 


T  =  T  -  ah 
o  p 


(2-33) 


Integration  of  Equation  2-32  from  mean  sea  level  to  the  altitude  in  question 
yields 


G 

fa 


-  / 


dT 


J 


or 


G 

Ra 


—  [  /_  I  =  In  I  /r 


(2-34) 


Equation  2-33  is  substituted  into  Equation  2-34  with  the  result, 


Jtn  1  —  |  =  —  Zn 

Ra 


T  -  ah 

3. _ E 


» 


I 
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or 


Differentiation  of  Equation  2-25  with  respect  to  P0  gives 


« 


» 


This  equation  can  be  written  in  the  form, 


T  R 
o 

P  G 
o 


gR 

P  G 
o 


h  . 
P 


(2-35) 


(2-36) 


(2-37) 


) 
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Substitution  of  this  result  into  Equation  2-24  and  use  of  AP  =  (Pmsi  -  P0) 
gives 


+ 


<P 


msl 


gR 

P  G 
o 


(2-38) 


If  the  constants  are  now  evaluated,  Equation  2-38  becomes 


H  =  h  WP  -  29.92  in  Hg)  [924.87  ft/in  Hg 
ph  p  msl  L 

-  .006358831/in  Hg  (h  )]  .  (2-39) 


Equation  2-39  now  agrees  with  Equation  2-23  except  for  a  small  difference 
in  the  seventh  decimal  place  in  the  third  constant. 


The  constants  appearing  in  Equation  2-39  are  not  grouped  in  a  manner  that 
makes  their  aggregate  value  independent  of  the  unit  system  chosen.  In 
Equation  2-38,  P0  can  be  factored  out  to  give 


+ 


(2-40) 


where 

T  R 

H  =  -2_  and 
G 


A  = 


aR 

G 


The  constant,  H,  (which  has  dimensions  of  length)  is  called  the  Scale  Height 
and  is  defined  to  be  the  height  at  which  the  pressure  of  an  isothermal 
atmosphere  at  temperature,  TQ,  decays  to  e"1  of  its  surface  value.  The 


I 


♦ 


* 
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non-dimensional  constant,  A,  represents  the  temperature  lapse  rate  made  non- 
dimensional  by  scale  height,  H,  and  standard  surface  temperature,  T0.  These 
results  are 


H  =  518.688 


1716.5  ft2 
sec2  °R 


sec 


32.17405  ft 


or 

H  =  27672.24  ft 

and 


or 

A  =  0.1902562. 

The  three  constants  appearing  in  Equation  2-40  are  dependent  only  on  the 
units  used  for  the  standard  atmosphere.  The  equation  itself  is  not  aircraft 
dependent:  its  form  will  change  only  if  the  standard  atmosphere  is  changed  or 
the  underlying  assumptions  used  in  the  derivation  are  changed.  The  values  and 
physical  meaning  of  the  constants  appearing  in  Equation  2-40  are  given  in 
Table  2.2. 


» 


» 


I 
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Table  2.2  LOCAL  CONSTANTS 


Physical  Meaning  Dependencies 

temperature  lapse  rate  made  standard  atmosphere 

non-dimensional  by  scale 
height  and  standard  surface 
temperature 


scale  height  of  the  standard  atmosphere 

atmosphere 


Value/Units 

0.1902562 

unitless 

27672.24  ft. 
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SECTION  3 

ANGLE  OF  ATTACK  (AOA) 


* 


t 


3.1  General 

The  Angle  of  Attack  (AOA)  calculations  are  performed  by  converting 
the  input  value  of  the  indicated  AOA  to  its  true  value  and  smoothing  these 
converted  values. 

3.2  AOA  Scaling  and  Smoothing 

Background 

In  Reference  1,  the  OFP  listing  equation  for  final  smoothed  AOA,  a, 
is  given  as  (read  as  meaning  "set  equal  to”) 

a  =  a  =  (a  -  0.07842  a)  +  +  aIN^  ’  (3-1) 

a 

where  a  is  a  smoothed  AOA.  It  is  given  that 

a’  =  0.03921[ (1.598155x10-s  )a  -(0.0241111)1 

P  1  P 

and 

a  '  =  0 . 0392 1 f (1 .598155x10-5  )a  -  (0.0241111)1, 

IN  1  IN  J 

where  af)  and  ot JN  are  the  two  different  Indicated  AOA  sample  values. 

As  can  he  seen,  Equation  3-1  is  rife  with  numerical  constants  and  confusing 
notat ion. 


i 


The  AOA  system  consists  of  a  movable  vane  located  on  the  left  side 
of  the  fuselage  near  the  nose  of  the  aircraft.  The  vane  senses  the  angle 
between  the  free  stream  movement  of  the  air  (the  velocity  vector)  and  an 
arbitrary  fuselage  reference  line.  The  AOA  vane  is  immersed  in  the  turbulent 
flow  about  the  aircraft  and,  having  a  much  lower  effective  mass,  responds  more 
quickly  to  inputs  than  the  aircraft.  The  quick  response  of  the  vane  produces 
noise  on  top  of  the  actual  AOA  signal  produced  by  the  transducer.  The  signal 
must  be  filtered  to  remove  the  noise,  and  it  must  be  scaled  and  offset  to 
indicate  the  appropriate  angle. 

Definition  of  Variables 


We  define  the  following  variables  and  general  constants: 

T  computation  time  period 

a  represents  or  otp 

a  final  smoothed  AOA 

a 

a  smoothed  AOA 

a i  two  input  AOA  values  (ajN  and  ap)  averaged  over  fixed 

sampling  period 

ajN  first  indicated  AOA,  scaled  and  formatted,  within  fixed 
sampling  period 

ap  second  and  final  indicated  AOA,  scaled  and  formatted,  within 

fixed  sampling  period 

aj  true  AOA 

a  IN 


true  AOA  times  a  constant  (see  a^) 


a1  true  AOA  times  a  constant  f see  a  ] 

P  P 


aTll  indicated  AOA,  formatted  fseeaTtll 
IN  v  IN' 


a  indicated  AOA,  formatted  f  see  a  ) 

P  P 


Derivation 


The  transfer  function  of  the  AOA  system  is  given  to  be 


TS  +  1 


(3-2) 


where  x  is  a  time  constant  that  has  been  empirically  determined.  Equation  3-2 
can  be  written  as 


as  =  t 


(«!  -  a) 


Since  the  inputs  to  the  AOA  system,  am  and  Op,  are  sampled  over  a  period, 
T,  application  of  a  Z-transform  is  suitable  (see  Reference  2,  Appendix  A;  also 
see  References  3  and  4). 


(3-4) 


where  z  is  actually  the  Z-transform  operator.  Clearly  we  have 


-1  -sT 
z  =  e 


(3-5) 
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Since  an  exponential  can  be  written  in  the  form, 


we  take 


£ 

j  > 0 


-i 

z 


sT  , 


or 


s  = 


(1 


Thus,  Equation  3-3  can  be  written 


Actually  applying  the  z-operator  in  Equation  3-9  yields 


S3  S3 


In  view  of  how  is  defined,  we  multiply  both  sides  by  (1 
Equation  3-11  can  then  be  written 


»  G 


or 


-1 


-1 


(3-13) 


) 


» 


t 


n 


a 


n 


(aIN 


To  get  Equation  3-13  to  include  the  arithmetic  for  indicated  AOA,  we  note  that 
the  OFP  is  supplied  AOA  from  the  AOA  vane.  The  geometry  of  the  AOA  vane  is 
illustrated  in  Figure  3-1. 


From  Figure  3-1,  we  see  that 

a (electrical  AOA)  =  a^+  0  . 


(3-14) 


The  position  error  equation  has  been  calculated  by  systems  engineers  to  be 


°T  =  ‘l  V 


(3-15) 


Substituting  the  value  of  indicated  AOA  from  Equation  3-14,  we  get 
a_  =  k.  (a  -  0 )  +  k  , 


(3-16) 


or 


«T  =  k^a  -  k^9  +  k2  . 


(3-17) 


?. 


We  need  Equation  3-17  to  give  AOA  to  the  Arbitrary  Datum  Line  (ADL)  and  not 
waterline  -  100.  Therefore,  we  need 


aT  =  kf  -  4  , 


(3-18) 


where 


k2  =  ki9  *  k2  ‘  * 


» 
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Figure  5-1.  Diagram  of  Angle  of  Attack  Vane 


The  Signal  Data  Converter  (SDC)  scales  and  formats  the  DC  voltage  input  by  an 
angular  constant,  X.  Further,  the  AOA  is  scaled  by  an  angular  constant,  w,  to 
effect  a  conversion  to  OFP  angular  coordinates.  Thus,  Equation  3-18  becomes 


k2 


from  which  we  have 


and 


“IN  =  S  “IN  •  K2 


K.  a 
1  P 


-  K. 


where 


and 


k^X 


k ' 
K2 


(3-19) 


(3-20) 


(3-21) 


With  the  results, 


and 


(3-22) 


(3-23) 
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Equation  3-13  may  be  written 
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♦ 
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a  = 


(a  -  <a)  +  (a 


IN 


+  ap) 


where 


The  values  of  the  constants  are  given  in  Table  3.1. 
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Table  3.1  LOCAL  CONSTANTS 


Symbol 

Physical  Meaning 

Dependencies 

Value/Units 

kl 

none 

AOA  sensor  and 
test  data 

0.785525,  0.76 
unitless 

k2 

none 

AOA  sensor  and 
test  data 

0.4  degrees 

T 

time  constant 

function  of  weapons 
rate 

0.47007  sec 

T 

weapons  rate 

function  of  filter 
constant 

0.04  sec 

0 

angle  between  WL-100  and 
zero  resistance  of  electrical 
moment 

AOA  vane 

8  degrees 

0 

angle  between  WL-100  and 
pot  wiper  arm 

AOA  vane 

9  degrees 

* 

angle  between  ADL  and 

WL-100 

AOA  vane 

3  degrees 

X 

scales  DC  voltage  input 

SDC 

30/4096  degrees 

w 

conversion  constant  to 
change  units  to  circles 

units  to  be 
converted  are 
in  degrees 

360  degrees 
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SECTION  4 
DOPPLER  RADAR  SET 

4.1  General 

The  Doppler  Radar  System  is  a  sel f-contained  dead-reckoning  system 
that  obtains  desired  navigational  information  involving  aircraft  velocity  by 
means  of  a  Doppler  Radar  and  direction  by  means  of  a  directional  sensor. 

4.2  Corrections  to  Doppler  Ground  Speed 

Background 


Page  DRS-5  of  Reference  2  on  the  Doppler  Radar  Set  gives  the 


equation, 


Pl:  =  ground  -  ( .00025*ground*alt ) / 1000. 


(4-1) 


This  equation  corrects  the  Doppler  ground  speed  to  account  for 
earth  curvature  and  other  empirically  determined  inaccuracies  which  are 
functions  of  altitude. 


Figure  4.1  illustrates  flight  over  curved  earth.  The  aircraft  is 
located  at  0,  with  lift,  drag,  and  weight  vectors  as  indicated.  The  flight 
path  makes  an  angle,  y,  with  the  horizontal,  and  the  aircraft  moves 


along  the  path,  OP,  The  line  joining  0  to  0e  (the  center  of  the  earth) 
sweeps  along  the  earth's  surface  with  a  velocity,  X  (the  ground  speed).  The 
measured  value  of  velocity  in  the  aircraft  frame  of  reference  is  greater  than 
that  described  at  the  earth's  surface.  Thus,  the  measured  value  of  the  ground 
speed  must  be  corrected  for  aircraft  altitude. 

The  Doppler  radar  system  contains  other  errors  that  are  functions  of 
aircraft  altitude  or  radar  slant  range.  The  aggregate  value  of  these  errors 
has  been  determined  by  flight  testing  the  aircraft  at  various  altitudes  and 
comparing  the  measured  velocities  with  photographically  derived  data.  The 
curve  that  was  fitted  through  these  data  points  is  described  by  Equation  4-1. 

Definition  of  Variables 


We  define  the  following  variables  and  general  constants: 

A  point  of  intersection  of  line  0e0  with  the  earth's  surface 
D  drag  vector 

h  height  of  0  above  A 

L  lift  vector 

0e  origin  of  the  earth  centered  coordinates 

0  origin  of  the  aircraft  frame  of  reference 

r0  radius  of  the  earth 

V  magnitude  of  the  aircraft  velocity  vector  in  the  earth  centered 
coordinate  system 
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W  aircraft  weight  vector 

X  ground  speed 

Y  angle  measured  from  the  horizontal  to  the  velocity  vector 

Derivation 


The  equation  of  motion  for  flight  in  a  great  circle  plane  is 
•  r 

X  =  V  — —  cos  y  , 

r  +  h 
o 


(4-2) 


where  V  is  the  velocity  vector  of  the  aircraft  in  the  earth  centered 
coordinate  system.  If  y  =  o,  which  implies  flight  along  a  path  parallel  to 
the  earth's  surface,  Equation  4-2  becomes 


(4-3) 


The  right  hand  side  of  Equation  4-3  can  be  expanded  with  the  binomial  series 
to  yield 


Since  the  quantity,  h/rQ,  is  very  smali,  terms  of  second  and  higher  powers 
may  safeiy  be  neglected,  thus  resulting  in 


If  the  value  for  rQ  is  substituted  into  Equation  4-4  in  the  same  way  as  in 
Equation  4-1,  the  result  is 


X  =  V 


KhV/1000  . 


(4-5) 


» 


t 


The  value  given  in  Table  4.1  for  K  is  about  1/5  of  that  used  in  Equation  4-1. 
The  difference  can  be  accounted  for  by  the  fact  that  the  value  of  K  represents 
not  only  that  correction  due  to  flight  over  a  curved  earth  but  also  a  number 
of  empirically  determined  correction  factors.  To  proceed  further  with  an 
analytic  derivation  of  the  given  value  of  K  would  involve  a  detailed  analysis 
of  the  specific  Doppler  radar  on  the  A-7  aircraft. 


Finally,  we  need  to  pose  either  Equation  4-1  or  Equation  4-5  in  a  form  which 
contains  constants  of  the  correct  variable  type.  The  best  form  should  be 


(4-6) 


with  r0  =  20,9 30,000  ft. 


> 


I 


> 


I 


» 
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Table  4.1  LOCAL  CONSTANTS 


ymbol 

Physical  Meaning 

Dependencies 

Value/Units 

K 

correction  due  to  flight 

set  of  flight 

4.8473x10‘b  i 

over  curved  earth  plus 

test  data  used 

altitude  correction  factors 
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APPENDIX  A 


LIST  OF  SYMBOLS 


G  deceleration  due  to  gravity  at  mean  sea  level 

32.17405  ft /sec2 
9.80665  m/sec2 

g(z)  acceleration  due  to  gravity  as  a  function  of  altitude 

hp  geopotential  altitude 

hp  ADC-sensed  barometric  altitude 

H  scale  height  of  the  atmosphere  27672*24  ft* 

Hp^  pressure  altitude  corrected  for  actual  pressure  at  sea  level 

Mn  true  Mach  number 

PQ  standard  pressure 

29,92  in,  Hg. 

101332.27  N/m2 
2116.2  lb/ft2 

Ps  static  pressure 

Pt  total  pressure 

q  dynamic  pressure,  P^  -  Ps 
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LIST  OF  SYMBOLS  (Continued) 


R  the  gas  constant  for  dry  air 

1716.5  ft2 /sec2 -°R 
287.04  m2 /sec2 -°K 

T  computation  time  period 

T0  standard  temperature 

518.688  °R 
288.16  °K 

Vvv  vertical  velocity  or  rate  of  change  in  altitude 

z  tape  line  altitude,  i.e.,  the  actual  distance  from 

mean  sea  level  to  the  point 

a  standard  temperature  lapse  rate  of  the  atmosphere 

0.00356616  °F/ft 
0.0065  °C/m 


a  final  smoothed  Angle  of  Attack 

smoothed  Angle  of  Attack 

a j  two  input  Angle  of  Attack  values  (ajN  anc*  <*p)  averaged 

over  fixed  sampling  period 
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LIST  OF  SYMBOLS  (Continued) 


t 


7 


* 


aTkl  indicated  Angle  of  Attack,  scaled  and  formatted 

IN 

dp  indicated  Angle  of  Attack,  scaled  and  formatted 

a j  true  Angle  of  Attack 

a '  true  Angle  of  Attack  times  a  constant 

IN 

a*  true  Angle  of  Attack  times  a  constant 

P 

oi  ,  indicated  Angle  of  Attack,  formatted 

IN 

a  indicated  Angle  of  Attack,  formatted 

P 

Y  ratio  of  specific  heats 

1.4  for  air 

t  time  constant  associated  with  the  static  pressure  system 

AP  difference  between  standard  pressure  and  that  existing 

at  sea  level 
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APPENDIX  B 


REAL  VARIABLE  TYPES 


Variables 

Units 

Acceleration 

1/ 12 

Angle 

unitless 

Angular  Rate 

1/t 

Density 

m/fc3 

Distance 

Z 

Mach 

unitless 

Pressure 

F/l2 

Speed  or  Velocity 

lit 

T  ime 

t 

I 
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APPENDIX  C 

Dependencies  of  Constants  used  in  Physical  and  Empirical  Equations 


t 


t 


t 


Suppose  that  the  "true"  physical  equation, 


y  -  f[ xi * • • • f*m;  ^ . kn) 


f or  1|  ^  K.  U|  (i-1 , , , .  ,m) 


(C-1) 


where  are  independent  variables  and  ki,,..,kn  are  physical 
constants  (the  actual  values  of  which  may  have  been  empirically  determined)  is 
known.  Then  each  k|  (i=1,.,,,n)  can  be  described  in  terms  of  its  particular 
physical  dependencies. 


On  the  other  hand,  suppose  Equation  (C-1)  can  be  numerically  represented,  over 
some  regions  (possibly  disjointed),  by  the  set  of  p  empirical  equations, 


=  ^1 jf • • • j) 

for  Ljj  _<  X|  _<  Ujj  (i=1 


(C-2) 


,m;  j=1,...,p) 


where  Kjj  ( i= j=1,..,,p)  are  empirical  constants  -  i.e., 


» 
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> 


I 


t 


» 


t 


I 


(C-3) 


( i  —  1 ,  •  »  •  ,m) . 

Then  for  any  r^  set  of  constants  K^r  ( 1=1 ,  *  • .  ,Nr) ,  each  and  every  one 
has  the  same  dependencies  as  any  other  -  viz.,  totally  dependent  on  the 
process  used  to  develop  the  equation 

Yr  =  Fr  ( xi > • *  *  >xm>  Kir> • • • »Knrr) 

(C-4) 

for  Ljr  _<  xi  _<  Uir  (1=1, ...  ,m). 

These  dependencies  would  include  characterist ics  of  measuring  devices  and  even 
the  particular  test  data  used  to  represent  Equation  (C-1). 

(It  is  worth  noting  that  if  we  consider  Equation  (C-2)  for  the  special  case  of 
p- 1  and 

Ni  =  n 
Kn  =  kL  <1=1, 

=  1 i  (1=1 , . . . ,m) 

uil  -  ui  U=1f  •••»«) 

that  this  resultant  "empirical"  equation  becomes  Identical  to  Equation  (C-1).) 


y  = 


Vx  for  Lu  <  <  Ui! 


V2  for  L|2  jC  Xi  _<  Ui2 


Vr  for  Lir  _<  Xi  _<  Uir 


Vp  for  Lip  _<  xi  _<  U 


IP 
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APPENDIX  D 


NOTATION 

While  most  readers  are  probably  familiar  with  the  mathematical 
notation  used  in  this  report,  what  follows  is  a  list  of  notational  symbols 


that,  due  to  inconsistent  usage  in  the 


x  =  y 

x  *  y 


2> 

R(j) 

0(f(x)) 


*For  a  detailed  discussion,  the  reader 
Computer  Programmi ng, "  Vol.  1,  2nd  ed. 


literature,  ought  to  be  clarified  here, 
x  is  identical  to  y 
x  is  approximately  equal  to  y 

identical  to  F(kj , • • • ,k  )  for 

K*i. ••••*„)  ■  f(*i> 

8xJ 

n  t 

binomial  coefficient  - : - 

j ! (n-j ) ! 

sum  of  all  f(j)  such  that  j  is  an 
an  integer  and  relation  R(j)  is 
true.* 

any  quantity  F(x)  such  that 
J  F  ( x )  |  <_  N  |  f  ( x )  |  over  the 
interval  a  £  x  £  6,  where  N  is 
an  unspecified  constant.* 


is  referred  to  Knuth,  D.E.,  "The  Art  of 
Addison-Wesley ,  Reading,  Mass.,  1973. 


